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1.   Introduction 

It  is  well  known  that  surface  waves  can  be  excited  and  propagated 
along  particular  idealized  electromagnetic  structures.   Examples  are  dielectric 
coated  surfaces^  dielectric  slabs  and  corrugated  or  laminated  structiures. 
In  so  far  as  the  surface  wave  feature  is  concerned,  it  is  known  that  the 
details  of  the  structure  ney  be  suppressed  and  replaced  by  an  impedance  boun- 
dary condition  [q]  .   Subsequent  developments  showed  that  the  use  of  impedance 
boundary  conditions  facilitated  the  solution  of  more  complex  problems  such 
as  the  excitation,,  propagation  and  diffraction  of  surface  waves  on  structures 
with  discontinuities  in  impedance  and  geometry  [5,6,  15-2Jj-]  .   The  pheno- 
menological  representation  by  an  impedance  boundary  condition  is  also  useful 
in  analyzing  the  effect  of  other  structures  such  as  absorbers  and  surfaces 
with  anisotropic  conductivity  [6,  l6] . 

Since  an  open  surface  structure  acts  to  some  e:ctent  like  a  wave- 
guide, and  since  its  theory  is  of  use  in  surface  wave  antennas,  it  would 
seem  to  be  of  interest  to  generalize  the  existing  theory  to  mult i- mode 
surface  wave  guides .   This  means  to  generalize  the  body  of  theory  mentioned 
above  to  the  case  of  surfaces  that  support  more  than  one  surface  wave.   Thus 
one  should  ultimately  possess  a  body  of  theory  for  surface  wave  guides 
analogous  to  that  already  available  for  conventional  wave  guides.   Typical 
problems  of  interest  are  the  analysis  of  the  effects  of 


*  The  Bibliography  given  at  the  end  of  the  paper  is  representative  but 
not  complete.   Additional  references  are  given  by  Collin  [7],  Cottony  [8] 
and  Zucker  [25] . 
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radlation  from  ,'    '    •       3  in  cross  section,  effects  -  ..s, 

etc.  'i  paper  is  the  first  part  of  a  series  .  ,    ;  to  these  topics. 

The  phenomenological  theory  on  which  the  generalization  is  based  has 
not  been  v   '  '    ;-cperimentally.   However,  it  is  clear  that  such  siorfaces 
can  be  realized,  at  least  by  sufficien  -.  trie  slabs  or  coatings. 

There  are  also  certain  qualitative  features  of  the  results  for  the  phenomenological 
theory  of  one  mode  in  the  case  of  the  excitation     '  :'m  which  are  in  agr 
with  the  exact  t:_     leal  results  for  a  slab  or  thin  coating.   It  would 
certainly  be  worthwhile  to  compare  the  exact  theory  with  the  phenomenological 

;  here  wit!      than  one  mode  and  in  detail.   Also,  exy      :tal 
results  for  one  or  more  modes  would  also  be  desirable.   Tlie  phenomenological 
theory  offers  the  possibility  of  analytical  solutions  of  complex  problems 
whose  analysis  from  the  standpoint  of  a  given  detailed  structure  is  difficult 
or  impossible  by  present  analytical  methods.  Another  problem  is  the  correspondence 
between  detailed  properties  of  the  structure  and  the  choice  of  a  generalized 
boundary  condition.   The  problem  has  been  treated  for  corrugated  structures 
which  support  one  surface  wave.  [7].   It  is  also  well  understood  in  the 
case  of  a  dielectric  structure  for  one  mode  and  has  been  partially  analyzed 
by  the  authors  (in  a  paper  to  appear  shortly)  for  a  multi-mode  dielectric 
structure.   In  brief,  it  would  seem  desirable  to  compare  the  phenomenological 
theory  (i.e,  exact  solutions  using  generalized  boundary  conditions)  with  exact 
solutions  in  which  the  detailed  structure  is  taken  into  account    ■      '   :ly. 
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The  purpose  of  the  present  paper  is  to  study  the  electromagnetic 
field  that  arises  when  a  magnetic  line  dipole  source  is  located  above  a 
plane  structure  that  supports  several  surface  waves  if  the  surface  of  the 
structure  is  characterized  by  a  generalized  impedance  boundary  condition.   This  con- 
dition is  introduced  here  for  the  first  time.   We  give  an  exact  solution 
for  plane  structures  supporting  one^  two  and  three  surface  waves.   We 
obtain  the  magnitudes  of  the  surface  waves  and  give  simple  formulas  for 
the  radiated  far  field  patterns.   We  also  show  how  the  methods  employed  in 
this  paper  can  be  extended  to  the  case  of  plane  structures  that  support 

any  number  of  surface  waves.  It  should  be  emphasized  that  the  effects  of 
terminations^,  finiteness,  etc.,  which  constitute  the  main  applications  of 

the  present  formulation,   are  excluded  from  consideration  in  this  introductory 

paper.   The  treatment  of  these  basic  diffraction  effects,  along  the  lines  of 

[5,  6,  15-2^] ,  are  the  subject  of  subsequent  papers. 

The  boundary  condition  we  employ  to  describe  the  propagation  of 

several  surface  waves  along  a  plane  surface  wave  structure  is  given  by 

n 

IT        (^  +     \.)u(x,y)    =  0      ,      y  =  0  (l.l) 

i=l         "^  ^ 

-H- 

where  u(x,y)   is  the   z-component   of  the  magnetic  vector   H(x,y),    x  and  y  are 

the  usual  cartesian  rectangular  coordinates,    and  ^.    are  constants  characteristic 

of  the    surface.      The  value   of  ^.    is   given  by 

X.    =  iue  z.    =  iwe(R.-iX.)    ,  (1.2) 

1  111  ^ 

where   e  is  the  permittivity  of  free   space,    to  is  the  angular   frequency,   and 

Z. ,    R.    and  X.    are  the   impedances,    resistances  and  reactances   of  the   surface, 
111  ' 

respectively.   For  surface  wave  propagation  it  is  necessary  to  require  that* 
Re  A.  >  0  .  (1.3) 


*  This  assumption  is  employed  in  the  illustrative  examples  that  follow, 
but  is  not  an  essential  restriction. 
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lt is  worth  noting  that  if  there  is  to  be  negligible  attenuation  of  the 
surface  waves  along  the  surface,  it  is  n      :-y  to  require  that  R^  <  <  X^. 
This  boundary  condition  is  a  natural  generalization  of  the  usual  impedance 
boundary  condition  employed  on  surfaces  supporting  one  surface  wave.   For 
example,  when  n  =  1  in  (l.l)  we  obtain  the  boundary  condition 

{^  +  \Mx,j)   =0  ,  y  =  0     .  (i.M 

This  boundary  condition  is  well  known  and  is  a  good  approxi;:-     -  to  several 
important  physical  configurations.   Perhaps  the  simplest  example  is  a  thin 
dielectric  layer  lying  on  a  perfectly  conducting  ground  plane.   This  and 
other  possible  physical  configurations  are  discussed  in  \j] ,   where  additional 
references  are  given.   It  is  clear  that  a  surface  wave  solution  of  the  form 

-^  y  +  i/k^+X^   |xl 
ujx,y)  =A^  e  -^         "^  (1.5) 

satisfies  the  wave  equation,  is  outgoing  at  infinity,  and  satisfies  (l.^). 
Hence  a  surface  characterized  by  (l.^)  can  indeed  support  sirrface  waves. 
Similarly,  the  generalized  boundary  condition  (l.l)  can  also  support      -  3 
waves  of  the  form 


o 

Ixl 


^—      -  ^ .  y  +  iy  k  +XT  |  x  | 
(x,y)  =  2^  A.e   ^         ^       .  (1.6) 


u 

1-1 


Therefore  (l.l)  is  the  natural  extension  to  str     ■  -s  supporting  several 
surface  waves.   We  shall  call  these  multi-mode  surface  wave  structures  and 
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use  the  generalized  boundary  condition  (l-l)  s-s  their  defining  property. 
We  will  not  discuss  the  physical  realizability  of  the  siirface  wave  structures 
here.   V/e  point  out,  however,  that  a  perfectly  conducting  ground  plane  with 
a  sufficiently  thick  dielectric  coating,  or  several  different  thin  dielectric 
coatings,  will  support  more  than  one  surface  wave. 

In  Section  2  we  study  the  electromagnetic  field  that  arises  when  a 
magnetic  line  dipole  source  is  located  above  a  plane  structui-e  that  supports 
a  single  surface  wave.   In  the  special  case  when  the  line  source  is  on  the 
surface,  the  m.ethod  and  result  reduce  to  those  given  previously  in  [20]  . 
Nevertheless,  the  treatment  is  complete  and  given  in  some  detail  since  the 
procedure  is  the  same  when  more  surface  wave  modes  are  considered.   The 
basic  ideas  for  obtaining  the  exact  solution  and  the  calculation  of  the 
surface  waves  are  contained  in  this  section. 

In  Sections  3  and  h   we  study  the  electromagnetic  fields  that  arise 
when  two  and  three  surface  vraves  are  excited  on  a  plane  structure, 
respectively.   The  fields  in  both  cases  are  produced  by  a  line  magnetic  dipole 
located  above  the  surface.   In  Section  5  "we  consider  the  case  of  any  number 
of  modes  and  summarize  the  main  steps  for  treating  this  problem.   In  addition 
to  the  surface  waves,  the  radiated  far  field  is  of  considerable  interest.   In 
Section  6  we  determine  the  far  field  radiation  patterns  for  one,  two  and  three 
modes  and  present  them  in  a  form  suitable  for  computation.   The  method  is 
direct  and  employs  an  elementary  procedure  that  does  not  require  the  exact 
results  obtained  in  the  earlier  sections.   We  also  show  how  to  find  the 
radiation  patterns  for  any  number  of  modes. 


2.   Single- Mode  Surface  V^ 


We  first  study  the  electromagnetic  field  tliat  arises  when  a  magnetic 
lin-  -'e  is  located  above  a  plane  structure  that  supports  a  single 

surface  wave.  Consider  the  plane  surface  y  =  0  and  suppose  a  magnetic  line 
dipole  is  located  at  the  point  x  =  0,  y  =  y^.  See  Figure  1.  In  the  region 
y  ^  0  V  we  have  free  space.   The  bound-  ,    edition  on  the  surface 

is  given  by  the  impedance  condition 

(^  +\)u(x,y)  =  0      ,   y  =  0  (2.1) 

where  u(x,y)  is  the  magnetic  field  component  H  and  ^  is  a  c    '     ■"   ■     '  ':ic 
of  the  surface  material.   We  wish  to  solve  the  time-re!;'     iixwell's  equations, 
subject  to  the  prescribed  condition,  and  obtain  the  amplitude  of  the  siirface 
waves  propagated  along  the  surface. 

The  time  d-      :t  form  of  Nb-xwell's  equations  is 


VxH  =  -iue  E 


Vx  E  =  itoji  ff 


(2.2) 


where  E  and  H  are  the  electric  and  magnetic  field  intensities,  and  e  and 

^  are  the  permittivity  and  magnetic  permeability  of  free  space.   We  assume 

-iiot 
the  time  dependence  to  be  of  the  form  e    .   Because  of  the      '  r-y,  the 

field  produced  is  independent  of  z  and  hence  is  ae 

value  of  H  =  u(x,y).   We  have 

H  =  H  =  E  =0  (2.3) 

X    y    z 


-6a. 


SOURCE 


/////////////////////  // 


MAGE 


ir.B) 


/////////////////// 


(o,-yo) 


►  X 


Figure  1.   IVfegnetic  Line  Dipole  Source  and  Image  in  the  Presence  of  a 
Surface  Characterized  by  a  Generalized  Impedance  Boundary  Condition. 
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(2.4) 


and 

1   ba 

_    1   lii 
The  function  u  =  u(x,y)  satisfies  the  wave  equation 

(^2+  k^)u  =  -lw5(x-0)5(y-y^)  (2.5) 

where  v  is  the  rectangular  Laplacian^  k  is  the  propagation  constant  of  free 
space  and  5  is  the  Dirac  delta  function.   Therefore,  the  mathematical  problem 
reduces  to  that  of  solving  the  inhomogeneous  wave  equation  (2.5)  subject  to 
the  mixed  boundary  condition  given  by  (2.1).   In  addition  to  the  prescribed 
boundary  condition,  we  require  the  far  field  to  be  outgoing,  and  the  value 
of  the  total  magnetic  field  (excluding  the  soirrce)  to  be  finite  everywhere. 

We  now  employ  a  technique  devel  ;     ;:  the  authors  ["5^6,16-22]]  for 
use  in  solving  mixed  boundary  value  problems  arising  in  electromagnetic 
diffraction  theory.   Let  us  first  introduce  an  auxiliary  function  v  defined  by 

V  =  (|p  +\)u  (2.6) 

Then  v  satisfies  the  equation 

(^  +  k^)v  =  -lirt(-  -^  +  \)S(x-0)5(y-yQ)  (2.7) 


o 


where  we  have  used  the  result  -r-^iy-'Y^)   -  -   ■^r-S(y-y^).   Also,  v  is  now  subject 
to  the  sim.pler  conditions 


(1)  V  =  0   ,    y  =  0 

(2)  Outgoing  traveling  waves  at  <x> 


> 


(2.8) 


If  the  problem  for  v  can  be  solved^  it  is  possible  to  find  a  solution  for  the 
problem  involving  u^  since  the  differential  equation  (2.6)  can  be  integrated. 
A  particular  solution  of  (2.6)  is 

Up(x,y)  =  e     J      e  v(x,Ti)dT]  (2.9) 


Tlie  solution  of  the  auxiliary  problem  is  easily  constructed  by 
introducing  the  image  source  located  at  the  point  x  =  0,  y  =  y    The 
desired  solution  for  v  is 


V 


where 


(x.y)  =  1^-  I;  +  \)     [4^h^^_)  -   H^^)  (kr^)}  (2.10) 


/^H^Xf    •  (2.11) 


Here  r  and  r  are  the  distances  from  the  source  and  image  to  the  observation 
point,  respectively.   Clearly  the  solution  vanishes  on  the  boundary,  is  outgoing 
at  infinity,  and  possesses  the  correct  singular  behavior  at  the  source  x  =  0, 
y  =  y  .   If  we  substitute  (2.10)  into  (2.9)  and  integrate  by  parts  we  obtain 


Up(x,y)=irt  |H^^^(kr_)  +H^^^(le-_^)1  -2niX^e  ^^^  e  ^  H^^^(M_^)dTi    (2.12) 
where  R  is  the  same  as  r  except  that  y  is  replaced  by  the  integration  variable 
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T).      The   funci.  ■  .(x,y)    just   found  satisfies  the  original  \,^--~~^-~-^-  -   --"tion 

(2.1),  3ses  the  correct    source  behavior  and  is   outgoing  at  00.      However, 

it   is  not  a  wave   function  since  the  x-deriArative   is  disco:,    '  :s  in  regions 

including  the  positive  y-axis.      This  defect  can  be  remedied,   as  will  be 
shown  below,   by  adding  to   (2.12)  another  '    n  of  the  wave  equation  as 

follows : 


u(x,y)      =  i«|H^^^(kr_)    +  H^^^(kr^)]       -   2^1^^^  ^^  f    e  ^"^  e  ^     H^^^(kR^)d 

—00 

-X^y    +   i/?7I^|x| 


(2.13) 


+  L,^       .3 


The  introduction  of  the  additional  term  with  the  unknown  constant  A  now  makes 
it  possible  to  satisfy  the  continuity  condition  involving  -^  .   We  note  that 
the  form  of  the  adn '         •  represents  a  surface  wave.   Hence,  when  we 
impose  the  continuity  ..        .it  on  ^  and  eAraluate  A  ,  we      :.  the 
amplitude  of  the  surface  wave  propagated  away  from  the  source.   This  is  the 
dominant  part  of  the  field  at  large  distances,  as  will  be  seen  later. 
We  now  proceed  to  determine  the  value  of  A  .   The  c    :..  Ity 
requirement  of  -^  across  the  line  x  =  0  can  be  expressed  as  the  following 
jump  condition: 

[I]   "  "0  [|(-0,  y)-|(x<0,  ,)]   =0  (2.U) 

From  (2.13)  we  obtain 


.10- 


^         ,      ^       f      r.^  n.  .  Ay    /^     \^ 


S  =     ^^  S: 


'"^  ^  ^  —00  + 


/-2 2  -X^y+i/l-^+  ^^ 

+  iA  /k  +  ^        sgn  X     e 


(2.15) 


where  we  have  used  the  identity 

•|;  H^^^kr)    -  -kF|^^(kr)   eos   0       •  (2.16) 

Substituting   (2,15)   into    (2.1i+)^   and  noting  that  the   first  terms   for  •^ 
yield  a   zero  contribution  to  the   jump,   we  have 


(2.17) 


+  2iA^  /k'~+  K     &  J 

The  limit  of  the  integral  on  the  right  can  be  evaluated  by  observing  that 
the  main  contribution  occurs  in  the  neighborhood  of  t)  =  -y^.   A  straight 
forward  calculation  yields 


_oo  + 


Using  this  result   in    (2.I7)  we   find 

A^  =--         r=^     e     ^  °        .  (2.19) 

/k   +  X^ 
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Thls  is  the  desired  amplitude  for  the   surface  wave.      The  exact   solution  of 
the  problem  is  given  by 


i(x,y)    =  irt       jn^^^l-a-J    +  H^^^(kr^)j      -    2niX^e'    ^^  f    e     ^"^   H^^^(kR^) 


dTi 
(2.20) 


+       ,  o  c 


Note  the  form  of  the  integral 


.X  -v      XX 


I  =\e  1  J     e^     H^^^  (k/x^+  (n+yj^hr]    .  (2.21) 

Tliis  integral  is  basic  in  the  present  paper.   In  fact,  the  solution  for  the 
multi-mode  problem  is  a  linear  combination  of  the  terms  in  (2.20)  with 
different  decay  constants  ^.  in  the  integral  and  surface  wave. 

3.   Two- Mode  Surface  Wave  Structure 
We  next  study  the  electromagnetic  field  that  arises  when  a  magnetic 
line  dipole  source  is  located  above  a  plane  structure  that  supports  two  surface 
waves.   Consider  the  plane  surface  y  =  0  and  suppose  a  magnetic  line  dipole 
is  located  at  the  point  x  =  0,  y  =  y^.   See  Figure  1.   The  boundary  condition 
on  the  surface  is  given  by  the  impedance  condition 


^¥^'i^^i"'-'2^'-'(^--y^  =0         ^   y  =  o  (3.1) 

where  u  is  the  magnetic  field  component  H^,  ^  and  >^  are  constants  characteristic 
of  the  surface  material^  and  ^   =f=  X    We  wish  to  solve  the  time-reduced 
Maxwell's  equations,  subject  to  the  prescribed  condition,  and  obtain  the 
amplitudes  of  the  svurface  waves  propagated  along  the  surface.   The  magnetic 
field  component  H^  =  u(x,y)  satisfies  the  wave  equation  given  by  (2.5). 
In  addition  to  the  prescribed  boundary  condition  (3.1),  we  require  the  far 
field  to  be  outgoing  and  the  value  of  the  total  magnetic- field  (excluding 
the  source)  to  be  finite  everywhere.   Note  that  once  u  is  known,  the  electric 
field  components  can  be  determined  from  {2.k). 

We  now  follow  the  same  procedure  as  was  used  in  Section  2, 
Let 

--(^-\)(^-\)--  (3.2) 

Then  v  satisfies  the  equation 


(V^+  k^)v  =  -hn 


— ^-(X+X)4-+^^ 


5(x-0)&(y-yQ)  (3.3) 


subject  to  the  simpler  conditions 

(1)  V  =  0  ,   y  =  0 

(2)  Outgoing  traveling  waves  at  =» 

A  particular  solution  of  (3.2)  is 


(3.M 
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Up(x,y) 


-X 


2   1 


^1^  a"  "-^ 

e    v(x,Ti)dTi  +  |— ^ 
l"  2 


e   v(x,Ti)dii. 


(3.5) 


The  solution  of  the  auxiliary  problem  is  easily  constructed,  as  before,  by 
introducing  the  linage  of  the  source.   The  desired  solution  for  v  is 


V 


(x,y)    =     irt 


r   ^ 


.:) 


_  (\  -t  xj  4—  +  ^  >^ 


H^/^(kr_)-   H^^^(kr^) 


(3.6) 


where  r       is   given  by   (2.11).      If  we   substitute    (3.6)   into    (3.5)^   and 
integrate  by  parts,   we   obtain 


Up(x,y)   =     irtf    H^^^krJ   -   H^^^(la^^)] 


5rtiX^(\^+XJ       _X^y 


e"^'^    H^^^(kR^.)dTi 


2ni^JX   +XJ      _X^y     /     X^n 


'X   -   Y 

1 


/I) 


e   '^   H^   ^(kR_^)dT 


(3.7) 


This  solution  is  not  a  regular  wave  function  since  the  x-derivative  is 
discontinuous  in  regions  including  the  positive  y-axis.   This  defect  can  be 
remedied  by  adding  to  (3.?)  another  solution  of  the  wave  equation  as  follows; 


-lil- 


u(x,y)   =  ±n  I  H^^^krJ   -   ni^\l^j] 


'O      '^^>j 


1  2  '^-oo 


1  2  ^-oo 


(3.8) 


-X^y  +  i/k'^+  X^    |x|  -^2^  ^  Vk'^+  ^2    '^-1 


+    A  e     ~  +     A  e 


The  introduction  of  the  additional  terms  with  unlmown  constants  A     and  A^ 
now  make   it  possible  to   satisfy  the  continuity  condition   (2.1J+).      From   (3-8) 
we   obtain 


ox 


I:   {h(^'(...,  -  H(^)(.n,,} 


2        1  ^-c»  + 


2rtikX    (X   +  \    )        _X   Y         X       Xn       .    , 
^  ,:x"        -        e      -         J        e   2      HJ^^kR^)      ^     d, 

12  ^-00  + 


(3.9) 


-X    y-H./k^+X^     |x|  ,^5 ^  -X    y-hi/k^+X^     |x| 


r2     2  ---^j^.^/^-  .-j_  K-i  /-o    72^  '-2-^ 

+  iA  /k   +  X        sgn  x  e  +  iA  /k   +  X        sgn  x  e 


Substituting    (3.9)    into    ('-.l^l),    and  using    (2.18),   we   obtain 
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^'«i\  (\+^J  -\yo 


and 


\   =      .,    ^ e  ^  "  (3.10) 


JmiA,(X  +X^)     .X  y 
Ag  =   .   -  -^  ^ e    °  (3.11) 


/k^+xl  (yxj 


These  are  the  desired  amplitudes  of  the  surface  waves.   The  exact  solution 
of  the  problem  is  given  by  (3.8),  where  A  and  A  are  given  by  (3. 10)  and 

(3.11). 

k.      Three-Mode   S'urface  Wave   Structm^e 

We  next   study  the  electromagnetic   field  that  arises  when  a 
magnetic   line  dipole   source  is  located  above  a  plane   structure  that   supports 
three   surface  waves.      Consider  the  plane   surface  y  =  0  and  suppose  a  magnetic 
line  dipole  is  located  at  the  point   x  =  0,   y  =  y   .      See   Firjio-e  1.      The  boundary 
condition  on  the   surface  is  given  by  the  impedance  condition 

where  u  is  the  magnetic   field  component   H     and  ^^,    ^     and  '^     are  constants 
characteristic   of  the   surface  material.      We  assiime  ^     /=  ^^  ;^  ^^.      We  wish 
to   solve  the  time-reduced  Nfexwell's  equations,    subject  to  the  pre, 
condition,   and  obtain  the  amplitudes   of  the   surface  waves  propti,_^^    .    ._      -  ng 
the   surface.      The  magnetic   field  component  u   satisfies  the  wave  equation   (2.5). 
In  addition  to  the  prescribed  boundary  condition   (-'kI),   we  require  the   far 
field  be  outgoing  at  infinity  and  the  value   of  the  total  magnetic   field 
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(excludlng  the  source)  be  finite  everywhere. 

Following  the  same  procedure  used  in  Sections  2  and  3?  we  let 


V  =  (4-  +  ^  )(|-  +  ^J(4-  +  ^  ) 


u. 


(i^.2) 


Then  v  satisfies  the  equation 


(V^  +  k  ")v   =  -hr 


a 


—   +  fX    +X    +X    )   ii—  -    (\  X^+X    X   +  X  X    )  ^ 
3  ^12      3^^.2        ^12        13        23^    ay 


^: 


+  X  X  X 
12   3 


^: 


6(x-0)5(y-yQ) 


0 
()+.3) 


subject  to  the  simpler  conditions 


(1)  V  =  0    ,     y  =  0 

(2)  Outgoing  traveling  waves  at  <» 


(i+.i^) 


A  particular  solution  of  (^-2)  is 


.X  V 
1-^ 


Up(x,y)   = 


(V\)^V^^ 


e    v(x,T])dTi  + 


_X  Y 


rx  -X  )(x  -x^j 

^  1   2-^^  3   2-' 


-X  y 


y     X 


e   v(x,T|)dTi 


(^.5) 


"*■   ("X  -X  )'(X  -X  ) 
^1  3'^^  2  3'^ 


e  -"  v(x^il)dT], 
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The   solution  of  the  auxiliary  problem  is 


V 


(x,y) 


irt 


_   ^     +  (X   +X+X    )  ^^  _    (X  X   +X  X   +X  X    )  4- 
.S         ^12      3''-,2        ^12     112-^''^, 


L    ^0 


+  X   X   X 
1  2  3 


(^r 


|H^l\kr_)   -    4^^(kr^)' 


C^.C) 


where  r   is  given  by  (2.11).   If  we  substitute  ('i.  6)  into  (^.^3)  and  integrate 
by  parts  we  obtain 


Up(x,y) 


in    [h^^^(}'^_)   -    H^^^(kr^)] 

2ni.\^(X^+XJ(X^+X^) 
^2     1^^   3      1^ 


\y 


2niX^(X^-fXJ(X^^+X^) 

— (X  _x  ){x  _x  j   ■ 

1     23      2 


•^oy 


^i"!    fi) 

e  1      H^l^(kR^)dTi 


r^        ^oT      fl) 


C'.T) 


2rtil    (X   +X    )(X,^+X    )        -X^y 


^.1 


3      ^      3 2_il      3      3         /        g    1      H;;"^(kR^)dTl 


rx  -X  )(x  -X  ) 

K   2_     3-"^    2      3^ 


^'  f  \.^... 


This  solution  is  not  a  wave  function  since  the  x-derivative  in  regions  including 
the  positive  y-axis  is  discontinuous.   We  can  remedy  this  defect,  however,  by 
adding  to  (h.^)   complementary  solutions  of  the  form 
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/"o      2~  /~2      ~  /~2      'c 

^^y+i  /k"+>^-L    |xl  -X^y  +  i/k  +x;    |x|  ->^  y+i  /k  +1. 


Then 

u(x,y)    =  Up(x,y)    +u^(x,y)  (i+.9) 

and  the  combined  solution   (4.9)   can  be  made  to   satisfy  the  continuity  condition 
on  -^  by  suitably  choosing  the  constants  A^   A     and  A   .      We  now  determine 

4^  from   (i^-.T),    (4.8)   and   (4.Q),   and  substitute  into    (2.14).      We   obtain 

tex    (X   +XJ(X   +X    )  -X   y 

A      =  ^      ^      " ^--i- e      ^  °  (4.10) 


v5?'+^(V\K^-\) 


4TCiX    (X   +X    )(X   +X    )  -X^y^ 

/k^+X^    (X^.X2)(X^-X^) 

4:tiX    (X   +XJ(X^+X    )  .X   y 

A^   =  3      ^     ^ 3  g      3   0  (^^^2) 


TK?    (X^-X    )(X^-X    ) 


3^13 

These  are  the  desired  amplitudes   of  the   surface  waves.      The   exact   solution  of 
the  problem  is  now  complete. 
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5.   Mi-ilti-mode  Surface  VJave  Structures 

The  soliition  of  problems  with  more  than  three  surface  wave  modes 
can  be  carried  out  using  the  same  procedure  described  in  Sections  2  through  '+. 
We  now  summarize  the  main  steps  for  treating  surface  wave  modes  of  any  order. 

The  mathematical  problem  to  be  solved  requires  a  solution  of  the 
wave  equation 

(^  +k^)  u(x,y)  =  .knb{x-0)b{j-y^)  (5.1) 

subject  to  the  boundary  condition 

fr   (^  +  ^^.hix,y)   =0    ,  y  =  0  .  (5.2) 

In  addition,  the  far  field  must  be  outgoing  and  the  total  magnetic  field 
(excluding  the  source)  must  be  finite.   We  first  make  the  substitution 

y  =     fr   (^  +  A.)  u  .  (p. 3) 

'  '    dy    J ' 

j=l 
Then  v  satisfies  the  equation 

(V2  +  k^)  V  =  ~hn   IT   (-  4-  +  ^  .,)5(x-0)5(y-y  )  (5.J|) 

j=l     ^-^0    ^ 

subject  to  the  simpler  boundary  conditions 


(1)  V  =  0     ,        y  =  0 

(2)  Outgoing  traveling  waves  at  <» 


(5.5) 
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The  solution  for  v  is  given  "by 

v(x.y)=i.ft  ^-h.^'^s^     {h^'^0-_)-^o'^^^-V)]  •        (5.6) 

A  particular  solution  of  the  original  field  u  can  be  found  by  solving  (5.3). 
It  is 

n    n      ^      -^±^    F       ^■i''^ 


^ —  ;r-l    ^  ,1   i 
~(i^. 


P 

i="l 


This  solution  is  not  a  wave  function  since  the  x-derivative  in  regions  including 
the  positive  y-axis  is  discontinuous.   We  can  remedy  this  defect,  however,  by 
adding  to  (5.7)  complementary  solutions  of  the  form 


i/k^+^r  |x|  . 
(x,y)  =  ^  A.  -^  '        ^  (5. 


-X  V  + 
i-' 
u 
c 

i=l 


Then  the  combined  solution  given  by 

u(x,y)  =  u  (x,y)  +  u  (x,y)  (5.9) 

P  o 

can  be  made  to  satisfy  the  continuity  condition  on  -^  .   This  is  equivalent  to 
imposing  the  jump  condition 


Lim   |^(x>0,y)  -  |H(^.<o,y)|   :=  0  .  (5. 10) 
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The  requirement  (5,10)  is  sufficient  to  determine  the  constants  A.,  which 
are  the  amplitudes  of  the  surface  waves.   The  exact  solution  for  u(x,y) 
will  contain  the  sura  or  difference  of  the  source  and  image  (depending  on 
whether  or  not  there  are  an  odd  or  even  niunber  of  surface  wave  modes, 
respectively),  a  linear  combination  of  integrals  of  the  type  (2.21)  with 
different  ?^.'s,  and  of  simi  of  surface  waves  as  given  by  (5.8). 

6 .   Determination  of  the  Radiated  Far  Field 
We  now  determine  the  radiated  far  field  arising  from  a  magnetic 
line  dipole  located  above  a  plane  surface  that  supports  surface  waves.   Tlie 
method  is  a  direct  and  elementary  procedure  developed  by  the  authors  \jj,6,1.6-22\ 
and  does  not  require  the  exact  results  obtained  in  the  earlier  sections. 

Consider  first  the  case  of  a  structiu"e  that  supports  a  single  surface 
wave.   The  same  physical  configuration  and  notation  used  earlier  applies. 
See  Figure  1.    From  (2.6)  and  (2.10)  we  have 


v=(4-+x)u  (6.1) 


'^ "  1' 


and 


V  =  in(-  |^-^  ^,)  {l4'^(l-J-H^^^(kr,)}   .  (6.2) 

Here  r     and  r     are  the  distances   from  the   soiurce  and  image  to  the   observation 

-      + 

point,  respectively,  and  are  given  by 
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/""2  2" 


+ 


'0^ 


(6.3) 


Let  r  be  the  distance  measui'ed  from  the   origin.      Then 


^J^ 


2     2 
r   =  yx   +y 


(6.ii) 


and  the  far  field  expresrjions   for  r       feecome 


i  ^  ^  t  Yq  sin  0   , 


(6.5) 


provid.ed  y^«  r.      VThen  kr        ic   large,    (6.2)  becomes 


V 


(r,0)~ 


kr 


^   sin(ky  sinO)  -  k  sin  0  cos(ky  sinO) 


i(kr-f) 


{6.6) 


where  we  have  used  the  asymptotic  expansion  of  the  Hankel  function  for  large 
argument.   We  e:cpect  the  radiated  far  field  of  u  to  have  a  similar  form, 
except  that  it  will  contain  surface  wave  terms  which  are  appreciable  at 
large  x  near  the  surface  y  =  0.   Hence,  for  fixed  0  and  large  r,  we  have 


0  /^  0  or  ir  , 


(6.7) 


where   f(0)   is  uniuaown.      Substituting   {6.6)  and   (6.7)   into    (6.1),   and 
using  the  result 
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i-   =  sin  0  ~   +  —     C03   0  4^    ~  sin  0   '^ 


^ 


br 


^ 


br     ' 


(6.6) 


.-1/2 


we  obtain  f(0)   correct  to  terms  of  order  r     '     .      The  expression  for  the   radiated 
far  field  then  becomes 


u(r,0);:.    yg 


pV    sin(ky  sinO)-k    sinO  cosCky^sinO)]         i(kr-|) 
FiksinO   +  a  "| 


(6.9) 


The  radiation  pattern  of  the  far  field  is  easily  determined  from  (6-9).   We 

X 

define  the  pattern   function  F(0,   p—   ,   ky^)   to  be 


Ho,   ^,   ky^)   = 


-i(kr-f) 

'Kr  ^        k 

ort 


u(r,0) 


(6.10) 


Hence 


F(0,   -    .    ^Jq) 


I    ~   sin(ky  sinO)-sinOcos(kypSinO)J' 


^(^f    -  cin^C^ 


(6.11) 


Consider  ne>ct  the  case  of  a   structure  that   :■ 
waves.      The  procedure  is  the   same  as  that   just   described. 
(3.6)  we  have 


:  two  surface 
From   (3.1)  and 


V 


u 


(6.12) 


and 


v(x,y)   =  iir 


^         (X  +  X  )  4-     +  K- 


Il'^'CkrJ-   H'^'Cki-^)"!    .(6.13) 
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The  radiated  far  field  for  v(rjO)  is 


^(r,0) 


A^X^   3in(kyQSinO)-(-\_^+)t,-,)k   sinOcosCky^sinO) 

2      2,  r  ^^^'^^^-f  ^ 


k  sin  0  sin(ky  sinO) 


(6.1M 


We   ezvpect  the  radiated  far   field  to  have  the   form    (6.7).       Substituting    (6.7) 

-l/2 
and   {6.1h)   into    (6.12)^   we   obtain  f(o)   correct   to  terms   of  order  r     '  ""    . 

Tlie   eicpression  for  the  radiated  far  field  then  becomes 


o        o 


u(r,0)  ^7g 


-        [(>.^X^-k'^sin"0)sin(kyQSin0)-(?>-^+>v,,)k   sinO  cos(kypSinO)]      i(kr-|) 

(6.15) 


Tik  sinO+X  1  Fik  sinO+A  1 


The  pattern  function  is  given  by 


F(0,    ^,   if   '    ""^O' 


-  X     X 


2  1  ?  I ' 

sin  O)sin(ky^sinO)- (~   +  ~)sinOcos(ky^sinO)    ' 


M'l         k 


'0 


[(~f+  sin^oj  [(~f+  sin^o] 


(6.16) 


The  calculation  of  the  radiated  far  field  for  a   structure  that   supports 
three   SLurface  waves  can  be  carried  out   in  a   similar  way.      Using   (1+.  1),   the 
asymptotic   form  of   (h.6),   and   (6.7)  we   obtain 
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Kkr.^) 


-k   sinO    (  (X^X,^+X^\   +X^X  J_k   sin  0  j     cos(kyQGinO) 


(6.17) 


The  pattern  function  is 


XX  X 


\2  2^-1 

(~)    +  sin  0 


(r^)'  +  sin^O 
k 


-1 


XXX 
1     £ 


XXX 
1 


?  o 


ff-    ^--^  if   ^  k-    ^    ^"^'"^  i       -i"(kyQ^i^®) 


(6.18) 


■^     ^    •      -^     3   +  ^  _2)   _    sin^o  C   sinO  coG(ky_sinO) 


^k     k     "*"  k     k     "^  k     k    ■^ 


'0 


For  structures  that  support  more  than  three  surface  waves,  the  procedure 
for  determining  the  radiated  far  field  is  the  same.   The  only  equations 
needed  are  (5.3),  the  asymptotic  far  field  form  of  (5.6),  and  (6.7). 
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